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We use the symmetry-constrained low-energy effective Hamiltonian of iron-based superconductors 
to study the Raman scattering in the normal state of underdoped iron-based superconductors. 

The incoming and scattered Raman photons couple directly to orbital fluctuations and indirectly 
to the spin fluctuations. We computed both couplings within the same low-energy model. The 
symmetry-constrained Hamiltonian yields the coupling between the orbital and spin fluctuations of 
only the same symmetry type. Attraction in the B 2 g symmetry channel was assumed for the system 
to develop the subleading instability towards the discrete in-plane rotational symmetry breaking, 
referred to as Ising nematic transition. We find that upon approaching this instability, the Raman 
spectral function develops a quasielastic peak as a function of energy transferred by photons to the 
crystal. We attribute this low-energy B 2 g scattering to the critical slowdown associated with the 
build up of nematic correlations. 

PACS numbers: 74.25.nd, 74.70.Xa 


I. INTRODUCTION 

The discovery of iron-based superconductors (FeSCs) 
opened new avenues in the research of strongly corre¬ 
lated systems jip— Despite the diversity in crystallographic 
structure and chemical composition, all the FeSCs share 
several generic trends. FeSCs are multiband and mul¬ 
tipocket materials. According to angle-resolved photoe¬ 
mission spectroscopy (ARPES) the Fermi surface (FS) 
contains two or three hole pockets at the center of the 
Brillouin zone, the F point, and two electron pockets cen¬ 
tered at (tt, tt), the M point in two-iron-atom unit-cell no¬ 
tations. The underdoped compounds undergo structural 
tetragonal to orthorhombic transition at the temperature 
Tg followed by or coincident with the spin density wave 
(SDW) transition at Tsdw- The superconductivity sets 

in when the magnetism is suppressed by dopingSrJS or 
1112 

pressur o 1 . 

The interplay between the magnetism and su¬ 
perconductivity is manifest in the weak-coupling 
renormalization-group analysis of competing instabili- 
tiesJ^^— The interaction amplitude in the spin-density- 
wave channel is renormalized in a way similar to that for 
the usual renormalization in the particle-particle chan¬ 
nel that normally leads to Cooper instability. Above the 
Fermi energy Ep the two channels affect each other. As a 
result, the interpocket pairing interaction is enhanced by 
the spin fluctuations, which were suggested to drive the 
unconventional superconductivity with the order pa¬ 
rameter changing sign between the electron and hole FSs. 
In this picture low (high) doping makes the magnetic 
(Cooper) instability a winner in a competition at energies 
below Ep. It follows that the proximity of the magnetic 
and superconducting phases on a phase diagram is not 
accidental. Hence, the understanding of magnetic and 
structural transitions is instrumental for the description 


of the superconductivity. 

Most commonly, at the magnetic transition the con¬ 
tinuous 0(3) symmetry and the discrete time-reversal 
symmetry are broken. In the FeSCs the spin alignment 
is magnetic along one direction and antiferromagnetic 
in the orthogonal direction. Such stripe magnetization 
lowers, therefore, the discrete rotational symmetry 
of the lattice down to 02 . The possibility of breaking 
the O 4 symmetry without breaking the spin 0(3) and 
time-reversal symmetry was studied in the context of the 
structural transition, and the corresponding transition 
was referred to as being Ising nematic.— In this picture, 
below Tg the spin-correlation length increases in one of 
the symmetry directions and decreases in the other, and 
the magnetisms sets in with little or no delay. 

The prevailing scenario of the structural transition is 
electronic. Specific to FeSCs is a rather high degree of 
ab anisotropy in the electronic properties. The resistivity 
anisotropy pb/pa in cobalt-doped BaFe 2 As 2 is reportedi^ 
to reach values as high as 2 for cobalt concentration 
X ~ 0.03, whereas the maximal orthorhombic distor¬ 
tion for the parent material, a: = 0, is only about 0.36%. 
Moreover, in stram-controlled samples the derivative of 
{pb — Pa)/(pb + Pa) with respect to the strain shows a di¬ 
vergence at the interpolated mean-field temperature T* 
= 116K for the parent compound. — The T* so obtained 
is only 22K lower than the actual transition tempera¬ 
ture, Tg= 138K. The relatively small difference Tg — T* 
is due to the lattice fluctuations being suppressed under 
the conditions of fixed strain. Likewise, the optical re¬ 
flectivity is nearly divergent at the nematic transitionj^ 
All these findings are indicative of the dominance of the 
electronic degrees of freedom in the nematic transition. 

It is, in general, hard to disentangle different electronic 
fluctuation channels breaking the same symmetry. At 
present the dominance of either charge or spin degrees 
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of freedom in driving the structural transition is not set¬ 
tled. There are two schools of thought as to the origin 
of electronic nematicity. — In the orbital nematicity sce¬ 
nario the difference in populations nxz — nyz of the dxz 
and dyz hon orbitals is believed to be the primary cause 
of the nematic transition. — In another scenario it is 
the spin that drives the nematic transition , Let mi ^2 
be the two staggered (antiferromagnetic) magnetizations 
on the even and odd iron sublattices, respectively. The 
nematic transition occurs when the two spin sublattices 
lock, (mi • m 2 ) ^ 0.— The two alternatives are the posi¬ 
tive and negative (mi • m 2 ), resulting in two orthogonal 
stripe magnetizations, = mi ± m 2 . These are 

the spin arrangements ferromagnetic in the X{Y) direc¬ 
tion and antiferromagnetic in the Y (X ) direction in an 
Fe-only lattice. The magnetic perspective is supported 
by the NMR data showing a low-T Curie-Weiss-like up¬ 
turn of a spin-lattice relaxation rate 1/TiT,— >2^ and by 
the scaling between the magnetic fluctuations and soft¬ 
ening of the elastic shear modulus at the structural tran¬ 
sition.— 

In this paper we do not attempt to resolve the above 
controversy, but rather explore the consequences of the 
nematic fluctuations as observed in recent Raman experi¬ 
ments— r— Even though the region of the phase diagram 
contained between Ts and Tsdw is either absent or quite 
tiny, the dynamical nematic fluctuations revealed by Ra¬ 
man spectroscopy kick in far into the paramagnetic phase 
up to room temperatures. The Raman spectroscopy is 
essentially a dynamic probe of electronic correlations of 
prescribed symmetry . The photon scattered inelasti- 
cally leaves some of its energy with the crystal. Selection 
rules fix the symmetry of the excitation, while the energy 
difference between the incoming and scattered photons, 
the so-called Raman shift, determines the energy of the 
electronic excitations. 


II. RAMAN RESPONSE IN THE FOUR-BAND 
MODEL 


A. Band structure model 

In this section we discuss the phenomenological four- 
band model based on the work of Cvetkovic and Vafek— 
In this model constructed using the method of invariants 
due to Luttinger— the interaction of electrons with light 
is easily obtained using the standard gauge-invariant 
minimal-coupling procedure^. Here we neglect the cou¬ 
pling between the different layers and consider the crys¬ 
tal structure to be quasi-two-dimensional, [see Fig.[T][a)]. 
Generically, in FeSCs each layer contains the iron atoms 
forming a simple square lattice with the basis unit vec¬ 
tors X and Y. The pnictogen or chalcogen atoms form 
the checkerboard with even and odd sublattices above 
and below the iron layer. Above the SDW transition the 
unit cell contains two iron atoms with the basis denoted 
by X and y. 


(a) 


(b) 




+As 
X 

oAs •Fe 






e i 







M 


FIG. 1: (color online) (a) The unit cell of a quasi-two- 
dimensional FeSC contains two iron atoms and two pnictogen 
atoms such as As (or a chalcogen atom such as Se). The atoms 
above and below the iron layer are denoted by crosses and by 
circles respectively. The basis vectors of the iron-only lattice 
are denoted by X and Y . The vectors x and y are chosen as a 
basis vectors of the two-iron unit cell lattice, (b) The two-iron 
Brillouin zone. The F point hosts two hole pockets and the 
M points hosts two electron pockets. The solid (black) and 
dashed (blue) lines denote the dxz and dy^ orbital contents 
respectively. The admixture of the dxY orbital at the outer 
parts of the crossed Fermi pockets at M is neglected. The po¬ 
larization vectors and e'^ for the B 2 g Raman configuration 
are shown. 


It is sufficient to consider a slightly simplified version 
of the model^ in which the four-dimensional effective 
Hamiltonian describing the M point is replaced with the 
two-dimensional one and the remaining electronic bands 
that are not crossing the Fermi level are discarded. We 
write, for the quadratic part of the Hamiltonian, 


n = 




k,Oi i,j —1,2 


+ fi^ka^kyi,j fj ,ka , (1) 


where cAc ifika) creates a hole (an electron) in a state 
with spin index a and momentum k counted from the F 
(M) point. The index i = 1 {i = 2) refers to the dxz 
(dyz) orbital content. For that reason the Hamiltonian 
o refers to the orbital basis and reads 


ni = 


er + 


2mr 

ikl-kl) 


dkx k' 


xnyy 


er 


m-ki) 
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( 2 ) 


for holes and 


o/M _ 
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€M • 


2mM 
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bkxky 
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( 3 ) 


for electrons. The parameters entering Eqs. © and 
([3]) obtained from the fits to the tight-binding calcula¬ 
tions^!^ are tabulated in Ref. ds^ . Below we set a = c, 
which corresponds to circular hole FSs. In this work we 
neglect the spin-orbit coupling, and at F, fc = 0, the two 
Bloch states are degenerate. Equation ([3|) neglects the 
admixture of dxy orbitals, and the parameter b is the 
pocket ellipticity. 

The Hamiltonians (HI) and (|3)) describe the band struc¬ 
ture shown in Fig.[T][b). The band structure obtained by 
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diagonalization of these Hamiltonians contains two hole 
pockets at F with orbital content alternating between 
dxz and dyz with tt periodicity and two electron pock¬ 
ets at M. The electron pockets cross, and their outer 
parts contain an admixture of the dxY orbital. Here we 
neglect such an admixture while preserving the overall 
symmetry of the Hamiltonian. 


B. Raman coupling 


Raman scattering is a two-photon process. Its am¬ 
plitude contains one part which is second order in the 
dipolar interaction and the first order in the coupling 
via the effective mass tensor. Assuming that the base 
frequency is detuned off the dipole transitions, it is cus¬ 
tomary to ignore the dipolar coupling. Under these cir¬ 
cumstances the inelastic Raman scattering cross section 
as a function of the Raman shift w is proportional to 
the imaginary part of the retarded Raman susceptibil¬ 
ity [K^{q,u})y'. We compute it from the corresponding 
Matsubara correlation function of the Raman vertices. 


K{q) = {ff)q , 


(4) 


where the vector q = {q,iuJm) includes the spatial wave 
vector q, and Matsubara frequency cOm, and we denote 

{AB)q = /J’ (irexp(fWmT)(Aq(T)B_q(0)). The experi¬ 
mental situation corresponds to <7 = 0 in Eq. 0 , and the 
Raman susceptibility k^{0,uj) is obtained from K{q) by 
setting <7 = 0 and performing the analytical continuation, 
iuJm —^ w. Below in writing the Matsubara frequency ujm, 
we omit the subscript m for brevity. 

The expression for the Raman vertices 


h,Oi i,j—1,2 


is fixed by the Hamiltonian, as formulated by Eqs. o, 
0 , and 0 , as well as the polarization vectors of incom¬ 
ing and scattered photons, and , 


Q2^r(M) 
Sr pipS ^ 

2^ ^ dkxdkx ■ 


( 6 ) 


In this work we focus on the i? 2 g Raman configu¬ 
ration such that polarization vectors of incoming and 
scattered photons are = (A -|- Y)lx/2 = x and 
e'® = (y — X)lx/2 = y, respectively (see Fig. [T]). The 
reason for this is twofold. First, the buildup of the low- 
energy B 2 g Raman intensity upon cooling is the domi¬ 
nant feature observed experimentally above Tg Sec¬ 
ond, both orbital and the nematic fluctuations have B 2 g 
symmetry. Indeed, Eq. ([0]) in combination with Eqs. 0 
and (l3|) gives 


r 

r = a 



= b 


I 

0 



(7) 


Equation © shows that photons in the H 23 Raman con¬ 
figuration couple directly to orbital fluctuations. 


III. EFFECTIVE ACTION AND RAMAN 
SUSCEPTIBILITY 


We compute the Raman susceptibility as given by 
Eq. (H]) with the Raman vertex specified by Eqs. © 
and ©. In doing so we follow closely the derivation of 
Ref. (jiffl . To compute the Raman susceptibility we add 
to the quantum action the source term. 


Sj — -f J—ujXg 


( 8 ) 


and the Raman susceptibility [Eq. ©] is obtained by a 
functional derivative of a free-energy functional, 


k{uj) = 


S'^T[J] 

SJu,SJ-u. 


(9) 


computed at J^i = J-ui = 0 . 

Here we focus on the spin interactions for definiteness 
and comment on the role of the orbital fluctuations. In 
the purely magnetic scenario of the nematic transition 
we write the interaction in the form 




~ 2 ^^ ^ ^ ' Si qSi -q, ( 10 ) 

q i—1,2 


where the spin operator is diagonal in orbital index i. 


.. = EE 




lllfk,ip 1 


( 11 ) 


k 1=1,2 


where (Jap are the Pauli matrices. The stan¬ 
dard Hubbard-Stratonovich transformation 

amounts to the decoupling of the interaction 
term [Eq. dlOp ] via the stripe magnetizations, 

oc Efe4+9,i(2)a‘^“/3/fc,i(2)/3- The integra¬ 
tion over fermion operators results in an effective action 
that closely resembles that of Ref. (0), 

J a' 


“I [|Sxy(0)|^ -I- |Sxy(w)|^] -I- [XALJu^Xvi^^) + C.C.] 


iq 
|2 I 1 ^ 


( 12 ) 


where we introduced the notation 

^xy(w) = ^ — Aq_^_,_QA^q _f2] . 

(13) 

For w = 0, Eq. m describes the classical contribution of 
the nematic fluctuations to the Ginzburg-Landau free en¬ 
ergy, while Sxy(w) describes the quantum nematic fluc¬ 
tuation driven by the external source at the same fre¬ 
quency. 

We now comment on Eq. (USD. First, we omitted the 
term oc [(A'’'^)^ -|- (A^)^]^ responsible for the renormal¬ 
ization of the spin susceptibility Xg that iscrucially im¬ 
portant for the nature of the magnetic and structural 
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phase transition. Here we are not concerned with either 
the feedback of nematic fluctuations on magnetism or 
mapping out the phase diagram, and for that reason we 
do not include this term in the action keeping the spin 
susceptibility unrenormalized. 

Second, the last term in Eq. m describes the coupling 
of the Raman vertex to the spin-nematic order parameter 
oc (A^)^ — (A^)^ via the triangular Aslamazov-Larkin- 
like vertex \al evaluated in the Appendix]^ We have 
shown this vertex can be approximated by a frequency- 
and momentum-independent function. Most relevant for 
the present analysis is the weak temperature dependence 
of Xal for temperatures not exceeding the mismatch be¬ 
tween the electron and hole Fermi surfaces. Because the 
latter can be a few tens of meV, the above statements 
hold for most of the relevant temperatures. 

We further perform the second Hubbard-Stratonovich 
transformation by introducing the nematic field (f). We 
introduce the static nematic field (f)o and the quantum 
time-dependent nematic field (j){T) = which cor¬ 

responding to the two terms quartic in A^^^^ in Eq. CH) 

. The resulting quadratic action reads 


^ ^ + <^0)1 Af |2 -f (x/ - ^o)| A^|2 

<? 

+ [{4’-uj + J_^)^xf(w) + C.C.] 

25 9 


(14) 


(a) 


X(q,uj + fl) 





(b) 


•9 



FIG. 2: Feynman graphs illustrating results m and (HU. 
The pair of thick blue wavy lines at the left and right ends 
of both graphs represents incoming and scattered photons. 
The black triangles represent the Aslamazov-Larkin triangu¬ 
lar vertex giving rise to the coupling constant Xal of light 
to the spin-nematic order parameter, (A^)^ — (A^)^, com¬ 
puted in the Appendix s The double wavy red lines denote 
the spin susceptibilities xilt ^)- While the graph in (a) is not 
specific to the XY geometry, the attraction in the nematic 
channel g > 0 makes it necessary to include the ladder di¬ 
agrams shown in (b), giving rise to the quasielastic peak in 
B 2 g geometry. 



where the dynamical spin-nematic susceptibility 

( 18 ) 


has been introduced. 

The action (II3 is quadratic, and the functional deriva¬ 
tive in Eq. ([9]) gives, for the Raman susceptibility, 


k{uj) = 3|AalP 


T(cc) 

l-3ffT(w)’ 


(19) 


The action (HI is quadratic with respect to the stripe 
order parameters, which thus can be integrated out ex¬ 
plicitly. This procedure results in the effective action 

exp(—J^) = J dA^dA^ exp(—S') 

^ Det(T>i^)Det(f>iL), (15) 


see Fig. [2] for diagrammatic representation. To compute 
T(w) we use the standard finite-temperature Matsub- 
ara summation technique over the discrete frequencies 
followed by an analytic continuation to the real axis, 
iujn —?> w -f fO, to obtain the retarded spin-nematic corre¬ 
lation function T(*a;„) —?> T^(w). 

To this end, we evaluate the bare susceptibility (ITSl) by 
converting the Matsubara sum over H into the complex 
integral 


= exp(-:S_kd 

I 25 g 


where we have defined matrices in the Fourier space 

TX± = ±{(l)-cj + XscJ-ui)M^±{(j)^ + XscJu>)M-, (16) 

with x±^ = Sg^^qJuA,uj 2 {Xqi^ ± 0o) and M± = 

Sq^.q 2 Suii,uj 2 ±i^- With the help of the standard formula 
we convert the determinant into the trace of the loga¬ 
rithm, DetTX = exp(Trln2?), and expand the resulting 
effective action up to second order in the nematic fields 
and the source strength. Such an expansion amounts to 
the mean-field approximation that can be justified in the 
large-limit. In Ref. (IT^ the thermodynamic properties 
of the same model were shown to be reasonably well cap¬ 
tured in this approximation for iV = 3, and we employ it 
here as well and write 

J^=^ + ^ - 3<^gT(0) - 3 !(<).,, + XalJu.^ T(w) , 
25 9 

(17) 


T(*a;„) = 


q 


dz z 

— coth —xi-iz + Wm, q)xi-‘iz, q). 

( 20 ) 


The integrand has two branch cuts at Im(z -1- iuJn) = 0 
and Im(z) = 0, where the product of two x functions 
has breaks of analyticity. As a result of the analytic 
continuation process, we get 


T(a;)=E 


dfl H 

— coth — 
27r 2T 


[x^{q, n + U})lmx^{q, H) -b Imx^(g, ^)x^iq, ^ - w)] 

( 21 ) 


To make further progress we use the standard expression 
for the spin-correlation functio n'^^d^ , 

x{q,^m) = 


^-^ + {Q-q? + \^m\/l 


( 22 ) 
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where the important scale is 


L - 2. 

Ts 


(23) 


c is a constant, 7 is the Landau damping coefficient, and 
Q = (tTjTt). Separating the real and imaginary compo¬ 
nents, one finds 


Imx^(q -h Q, n) = x" 




_[Lx_ 

(1 -I- ’ 


(24) 


Itex«(9 + Q. n) = x' = eg" ^ ■ (25) 

The symmetry between the two stripelike spin-ordering 
arrangements is broken at the Ising-nematic-type transi¬ 
tion. Here we assume the mean-field critical exponent 
v = 1 / 2 , i.e., 

ar) - WTn/{T-Tn), (26) 


with Tjv being the mean-field SDW transition temper¬ 
ature and I being a microscopic length scale. We em¬ 
phasize that the mean-field transition temperature can 
be substantially lower than the observed SDW transi¬ 
tion temperature, Tn < Tsdw- Equation (| 2 ^ also 
shows that the critical behavior of the static suscepti¬ 
bility x(^m = 0, Q) = oc (T —Tvr)"^ is as prescribed 
by the mean field. 

We proceed by substituting (l24l) and (1251) into the gen¬ 
eral relation m to get, for the imaginary part of the 
susceptibility. 


„ [°° dx f°° dy ycoth{y/2t) 

y + w y — w 

_ (1 -b x)'^ + (y + w)^ (1 -b x)'^ + [y — 


(27) 


Here we introduced dimensionless variables x = (g^)^, 
y = Hts, t = Tts, and w = ujts- The above expression 
is general. In a view of Raman experiments, below we 
consider in detail the limiting case of high temperatures, 
Tts ^ 1 , that corresponds to the regime of essentially 
classical fluctuations. 

In the classical region when t ^ I, assuming not 
too high frequencies, T > w, one can approximate 
coth(j//2t) Ri ‘2-t/y. The double integral in Eq. (l27l) can 
be then evaluated analytically, 


T"(a;) = 


4 T 7 


In 


w 



(28) 


Since the integral decays at a scale ^ 4rs, the approx¬ 
imation made should be reasonable for all frequencies. 
Similarly, we evaluate the real part of the susceptibility, 

X _ I __ (29) 

(I + xY + y2 (1 + xY + {y + wY' ^ ’ 


The integral in Eq. (|2^ is logarithmically divergent at ul¬ 
traviolet. We therefore isolate the divergence in Eq. (Ei) 
by focusing first on a static limit a; = 0. Then the differ¬ 
ence is convergent and can be easily evaluated as 

T'(w) = T'(0)-^ arctan((jjrs/2). (30) 

IOtt 

To evaluate the static susceptibility T'(0) we split the 
y-integration range in Eq. (ESI) into two regions, y < t 
and t < y < Axg, where the scale A is the ultraviolet 
cutoff. Making approximations coth(j//2t) r; 2t/j/ and 
coth(?// 2 t) r; sgn(?/) in the two respective intervals, the 
resulting integrals can be easily evaluated with the result 

tY a 

Tf''(0) « arctan(TTs) -b In -, (31) 

where the arctangent can be further safely approximated 
by tt/2. With these results at hand we find from Eq. (fT9l) 
for the imaginary part of the Raman susceptibility 


k"(w) 


0 x2 __ 

i?^(tn, r)-b (47r^g7)^/^(n7, r) 


(32) 


Here we have introduced dimensionless frequency vj = 
uj/Tn and temperature r = T/Tn, and also two dimen¬ 
sionless functions. 


I{va, r) = — In 
w 


R{xu,t) = 1- g 


I-b 


2tt^(t — I) 


(33) 


7r7 


^ + 7 L- -r arctan . 

r-1 ^ 4 V27r7(T-l) 


zu 


(34) 


where renormalized coupling g = 3gT]s[P/A tt, decay rate 
7 = and cutoff L = ln(A/T) should be used 

as fitting parameters; see Fig. [3] for results. The most 
prominent feature of our results is the critical enhance¬ 
ment of the Raman susceptibility upon approaching the 
structural transition with the characteristic buildup of 
the quasielastic scattering. 


IV. CONCLUSIONS 

In this paper we investigated theoretically the low- 
energy Raman scattering in underdoped FeSCs. The 
gross feature of the data is the quasielastic peak that 
gains in intensity and softens down at cooling above the 
structural transition. The phenomenon is observed exclu¬ 
sively in B^g Raman geometry. The Lorentzian-like fre¬ 
quency dependence of the B 2 g Raman susceptibility de¬ 
scribes the relaxation dynamics with the relaxation rate 
given by the position of the maximum. The temperature 
dependence of the susceptibility indicates the freezing of 
the electronic B 2 g fluctuations at cooling. Such behavior 
is naturally associated with the tendency to long-range 
order which breaks the B 2 g symmetry. That is, the broad 
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FIG. 3: (color online) The modeling of the quasielastic peak 
in the Raman response function in accordance with Eq. dsn, 
where k," is plotted in units of for the following choice 

of fitting parameters from the bottom curve to the top one: 
T = 2 (red dashed line), r = 1.5 (black thin solid line), and 
T = 1.25 (blue solid thick line); g = 0.2, 7 = 5, L = 10. The 
peak grows when the structural transition is approached upon 
cooling. 

relaxation-like feature can be attributed to the critical 
slowdown associated with the approach to the discrete 
symmetry-breaking transition. Upon cooling, the system 
experiences locking in one of the two degenerate configu¬ 
rations related to the C 4 rotation for increasingly longer 
time intervals. 

To understand the origin of the quasielastic peak as 
it appears in Eq. (HU), note first that in the static 
limit the real part of the Raman susceptibility scales 
as ^ (T — 9)~^. The temperature scale 9 < Tg, which 
can be explained in terms of the coupling of the elec¬ 
tron nematic fluctuations and the orthorhombic lattice 
vibrations studied recently in Ref. These fluctua¬ 
tions add to the static nematic coupling constant, gst = 
g + 7 ^/( 7 °Here 7 is a nemato-elastic coupling 
constant and (7° is the bare value of the orthorhombic 
elastic constant. The static coupling gst determines Tg. 
Crucially, however, it is the dynamic rather than static 
nematic coupling constant g < gst that defines 9 since 
the lattice response function has different static and dy¬ 
namic limits^. Correspondingly, the difference Tg — 9 is 
expected to correlate with the reduction of Tg in strain- 
controlled samplesi^. This is indeed reported to be the 
case in the recent measurements^. 

Distinguishing between different contributions to ne¬ 
matic correlations remains a challenge. Nevertheless, we 
can deduce the low-energy scattering by making the rea¬ 
sonable assumption on the imaginary part of the bare 
response. Assume that at low frequencies it scales as 
^ w/r, with r being noncritical atT = 9 and hence being 


a weakly temperature-dependent relaxation rate. Then it 
follows from the denominator structure of Eq. m that 
at low frequencies k"{uj) ^ {T — 9 + We thus 

see that the relaxation rate is suppressed by a factor of 
T — 9 compared to the bare rate T. We conclude that 
the quasielastic scattering is the case of critical slowing 
down. 

The intraband processes alone are insufficient to de¬ 
scribe the large frequency width of the quasielastic scat¬ 
tering. Indeed, at zero momentum such transitions are 
forbidden, and the quasielastic peak is absent i'^^d^ The 
small-momentum intraband transitions restricted to ei¬ 
ther r or M points are gapped and cannot account 
for quasielastic scattering either. The excitation of two 
electron-hole pairs at momentum close to the antifer¬ 
romagnetic wave vector enables the relaxation of zero- 
momentum excitations by lifting the kinematical con¬ 
straints. It was argued, however, that the phase-space 
limitations make the contributions of such processes to 
the relaxation rate scale as a cubic power of the fre¬ 
quency difference of scattered and incoming photons. — 
As our calculations demonstrate, this suppression is rel¬ 
evant only at very low frequencies, and for relevant tem¬ 
peratures and frequencies the scaling is essentially linear. 
In that regard this situation is very similar to that in 
cuprates 4^ Even though cuprates are single-band rather 
than multiband materials, the processes that matter the 
most are confined to the vicinity of hot spots, i.e., the 
points on a FS connected by the antiferromagnetic wave 
vector. 

Note added. Recently, we became aware of 
arXiv:1503.07646 in which spin-nematic susceptibility 
was computed in a similar model. 
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Appendix A: Derivation of the effective action in 
terms of the stripe magnetizations 
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The calculations in this appendix are an extension of the corresponding derivations in Re. dil). We introduce the 
eight-component spinor 


= {c\. 


fef 


ft ft 


(Al) 


in the direct product of orbital and spin spaces. Upon the introduction of the stripe magnetizations, via the 

Hubbard-Stratonovich transformation the action takes the form 


= - + - I |A^(g)r + |A'^(g)| 


(A2) 


where we denote k = (k,ie), and the Green’s function is 


where 


The free Green’s function in Eq. 
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where the two hole and electron Green’s functions are two-dimensional matrices expressed in a standard way, 

-1 


gif = {ze-ni:^^^ +Ef 


T,M 


(A6) 


through the hole and electron Hamiltonians ([2]) and ([3]) with energies counted relative to the Fermi level. We further 
have 


- 

^k,k' — 


- / <5, 


k-\-qf 


0 0 A-^(g) 0 

0 0 0 A^(q) 

A'^(q) 0 0 0 

0 A^(g) 0 0 


The source term according to Eqs. and (|5]) reads, 

k k' ~ '^=ta;^fc,fc'^e±u;,e 


0 


(A7) 


(A8) 


where the two-dimensional matrices are defined by Eq. ([7l). In the presence of the source (lASp . it is necessary 
to keep the term of the third order in V in the expansion of the free energy, 

5[A^, A^, J^, J_^] = \Tr{g^Vf + ^Tr{goVf + ^Trig^V)^ + [A^(q)|2 + \A^{qf , (A9) 

which is of prime interest for us in this work. While even-order terms in Eq. (IA9I) were considered in detail in Ref. (flH l , 
here we focus on the third-order term of Eq. (IA9I) . We specifically determine the contribution to the effective action 
that is linear in the sources J±uj and quadratic in A^’^. Such terms have the form 


TriV^^-goV^goV^go) = aJ±^ 

+ bfj±^TT 
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'1 0 ■ 
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In obtaining Eq. (jAlOl) we used the explicit form [Eq. ([7])] of the Raman vertices and the block-diagonal structure 
of all three matrices given by Eqs. (TASI) . (IAtI) . and (IA8I) . Furthermore the trace operation over the spin indices results 
in a summation over the Cartesian coordinates of the vectors of the stripe magnetizations A^A labeled by the index 
























£ = x,y, z. The trace operation in Eq. (lAlOl) includes, in addition to the usual trace of the two-dimensional matrices, 
the summation over fermion frequencies and momenta (ie, k) and boson frequencies and momenta {iQ, q). 

The integral of Eq. (lAlOll over the fermion frequencies and momenta is convergent at the upper limit. As a result, 
assuming that the Eermi energy is much larger than the typical energy carried by the magnetic fluctuation, we can 
neglect the boson momenta and frequency, q and fl, in the arguments of fermion Green’s functions. Furthermore, since 
the energy splitting of electronic bands is normally smaller than the Fermi energy, we can neglect it and approximate 
the hole Hamiltonian by a scalar function. In terms of the Green’s function we have 


’k,€ 




-1 


(All) 


where is the energy of the hole band relative to the Fermi energy with splitting neglected. Similarly, by neglecting 
the ellipticity-related energy that is small on the scale of the Fermi energy we arrive at the scalar Green’s function 
for electrons. 


-•M 




-1 


(A12) 


where is the energy of the electron bands relative to the Fermi energy with ellipticity neglected. Glearly, the above 
approximations make the corresponding Green’s function denoted by ^ and scalar and allow us to rewrite 
Eq. (lAllI) as 


/ k ^ 


TTiV^^-GoV^GoV^Go) ^ J±^Exy{±u;) I a (GE,f C + ^ (& 


(A13) 


Identifying the last term of Eq. (na) with Eq. (IA13|) . we finally arrive at the expression for the triangular vertex. 


Aal = 


Jk L 


^fc,£ 




k,€ 


(A14) 


where the expression for Sxy{xi) is given in Eq. CSl). For definiteness, we evaluate the first term oc a in expression 
(IA14I) . For simplicity we assume = —^e + 5he = C Ih® density of states to be a constant vq for both 

electrons and holes. Under these assumptions the substitution of the Green’s functions in Eqs. (lAIIl) and (IA12I1 in 
Eq. (IA14p . followed by the integration over yields 


{Qle) C = E 


e„>0 L 
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‘liCji -\- dell 
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As the Matsubara frequencies are of the form e„ = 2x{n -\- 1/2)T, we obtain from Eq. (IA15I) 




, / ; 1 . deh \ 

V’ ( 2 -*' 


'27rTy 
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(A16) 


where 4’’{x) is the derivative of the digamma function. An expression similar to Eq. (IAI 6 II holds for the contribution 
of the coupling of light to electrons, and we obtain 


Aal = (a + b)^lva 

TTi 


^ ' 2 2ttT 


(A17) 


We note that the contributions of electrons and holes are additive. As deh is typically a few tens of meV, we expect 
the inequality T < deh to hold. As the asymptotic expansion lm[4'{l/2 — ix)] ~ Ijx holds already for x > 0.25, we 
can write for the Aslamazov-Larkin vertex 


Aal ~ (a -b 5) 


deh ' 


(A18) 


We conclude that Aal is insensitive to the temperature variation in the relevant temperature range and is suppressed 
with increasing mismatch between the hole and electron Fermi surfaces. This result is in agreement with the alternative 
calculation for a different model.— _ 
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